Absence of BCS Condensation of Atoms in a Repulsive Fermi Gas of Atoms 
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Possible Bose-Einstein condensation of molecules and Bardeen-Cooper-Schrieffer (BCS) conden- 
sation of Cooper pairs of atoms in an ultracold Fermi gas of atoms with a repulsive two-body 
interaction are studied by using the path integral representation of the grand partition function. 
From the self-consistent equations obtained in the present work for the order parameters, we have 
found that BCS condensation of atoms can't occur in such a Fermi gas of atoms and that the conden- 
sate observed experimentally is composed of condensed molecules or possibly of preformed Cooper 
pairs. To substantiate our conclusions from the self-consistent equations for the order parameters, 
the effective atom-atom interaction mediated through molecules has also been computed with the 
Foldy-Wouthuyscn transformation and has been found to be repulsive, which implies a net repulsive 
two-body interaction between atoms and hence the absence of BCS condensation of atoms. 
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The experimental realization [lj of an ultracold Fermi 
gas of atoms marks the beginning of a new era of the 
scientific research on degenerate quantum Fermi gases. 
Since then, Bardeen-Cooper-Schrieffer (BCS) condensa- 
tion of fcrmionic atoms and Bose-Einstein condensation 
(BEC) of diatomic molecules of these atoms as well as the 
crossover between the BCS and BEC regimes have been 
under intensive study both theoretically 513,0 ,0,0,13,0 
and experimentally 0, El El El El 07 with the 
BCS-BEC crossover regime being e asily reached with the 
magnetic- field Feshbach resonance [lOj. 

In the BCS-BEC crossover regime, the nature of con- 
densation, i.e., whether the observed condensation is 
BCS condensation of Cooper pairs that are bound-states 
of fermionic atoms in momentum space or it is Bose- 
Einstein condensation of molecules that are bound-states 
in real space or it is condensation of preformed Cooper 
pairs of atoms, is still controversial 5|. In the experi- 
mental study of the crossover regime, both the sign and 
the strength of the two-body interaction between atoms 
can be tuned with the magnetic-field Feshbach resonance. 
It was reported that BCS condensation of atoms was ob- 
served on both the BCS (attractive two-body interaction) 
and BEC (positive two-body interaction) sides Ell- 
in this Letter, we concentrate on the BEC side of the 
crossover regime and investigate the nature of condensa- 
tion on this side. The Fermi gas of atoms is described by 
the generalized Hamiltonian 2] 

H = 2(i/-^)^a£a p + ^(e fe -/x)4 CT c feCT 
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where 2(v — fi) is the offset energy of a molecule rela- 
tive to the chemical potential /x, 



ajj and a T 



are creation 



and annihilation operators of a molecule of momentum 
p and energy 2(i> — /i), c k and c\~ a are creation and an- 
nihilation opperators of a fermionic atom of momentum 
k, spin <t, and energy = h 2 k 2 /2m, U (> 0) is the 
strength of the two-body interaction between atoms, g is 
the effective atom- molecule coupling constant, and N is 
the average number of atoms. Here the spin up and spin 
down actually represent two different hyperfine levels of 
atoms. Here we have assumed that there are only two 
components in the gas. 

Notice that the number of atoms N in the system ap- 
pears in the prefactors of the two-body interaction and 
atom-molecule coupling terms due to the Fourier trans- 
formation of the operators into momentum space. The 
presence of N in these prefactors in conjunction with the 
proper change of variables to be performed below will 
give the self-consistent equations for the order parame- 
ters the proper system size dependence. 

The aim of the present work is to study the nature of 
condensation on the BEC side of the crossover regime 
with the Hamiltonian given in Eq. ifTjl. The nature of 
condensation is described by the corresponding order pa- 
rameter. Therefore, to study the nature of condensation, 
we first set up the self-consistent equations for the order 
parameters of all possible types of condensation. Here 
we are concerned with two types of condensation, Bose- 
Einstein condensation of molecules and BCS condensa- 
tion of atoms. The order parameter for Bose-Einstein 
condensation of molecules is the number of molecules in 
the lowest-energy level, and the order parameter for BCS 
condensation of atoms is the energy gap in the excitation 
spectrum of quasiparticles. Our criterion is that a certain 
type of condensation is said to occur if the corresponding 
order parameter has a nontrivial solution below a certain 
temperature. 

Our starting point is the path integral representa- 
tion 0,0] of the grand partition function Z = Tr e~@ H . 
The order parameter for Bose-Einstein condensation is 
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introduced in the path integral representation of the 
grand partition function, whereas the order parameter 
for BCS condensation of atoms is introduced through 
the Hubbard-Stratonovich transformation 

[H0]1. The 

path-integral representation of the grand partition func- 
tion corresponding to the Hamiltonian in Eq. (QJ is given 

by mm 



Z = 



D[ip*,Lp;%l)*,ip] T exp<- / drV ip*(r) 
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where c.c. denotes the complex conjugate of the proceed- 
ing term, ^*(t) and <y5 p (r) are complex variables corre- 
sponding to the creation and annihilation operators aj,(r) 
and a p (r) of a bosonic molecule, and V , fc CT ( T ) anc ^ V'fco-('r) 
are Grassman variables corresponding to the creation and 
annihilation operators ct (t) and Cfe ff (r) of a fermionic 
atom. The integration measure D[tp* , Lp, ip* ,ip] T is given 
by 
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with At = (3/M. 

The path- integral representation of the grand partition 
function in Eq. Q is in imaginary time space and it can 
be transformed into imaginary frequency space through 
Fourier transformations of the complex and Grassman 
variables 

t Ppi T ) = y^ e ~ WmT yp(^m), (4a) 

iuj-m 

^fea(r) - e~ iu>nT i>k„{iu n ), (4b) 

where iu> m — i2mir / (3 and iuj n — i(2n + l)iv//3, with 
m and n being integers, are Matsubara imaginary fre- 
quencies for bosons and fermions, respectively. After the 
Fourier transformations, the complex and Grassman vari- 
ables are denoted succinctly by ip p and ipka , respectively, 
with short notations p = (p,iw m ) an d k = (k,iu> n ) in- 
troduced. The BCS order parameter A p is introduced 
through the Hubbard-Stratonovich transformation per- 
formed with respect to the two-body interaction term 



that, for convenience of the performance of this transfor- 
mation, is recast into the following form 
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For the repulsive two-body interaction term, the 
Hubbard-Stratonovich transformation must be imple- 
mented by making use of the Gaussian integral of the 
form e -<^/ Q = (a/2ni) j dz*dz e -«***+C**-C** with z 
and £ complex variables. That the two decoupled terms 
have different signs is because the two-body term has a 
minus sign on the exponential in Eq. (J2J . If an attractive 
two-body interaction term were to be decoupled, the two 
decoupled terms would have identical signs. Here the 
auxiliary field introduced in the Hubbard-Stratonovich 
transformation is denoted by A p which is the order pa- 
rameter for BCS condensation of atoms. It has been 
found that, to derive self-consistent equations for ip p and 
A p with a proper system size dependence, it is necessary 
to make a change of integration variables ip p — > ^/N<p p 



and A p -> \/N(3UA p . 

For the purpose of studying Bose-Einstein condensa- 
tion of molecules and BCS condensation of atoms, it 
is sufficient to retain only the zero momentum and fre- 
quency terms in the atom-molecule coupling term and in 
the decoupled terms arising from the two-body interac- 
tion. Then, the integration over the order parameters ip p 
and Ap with p ^ can be performed and the summation 
over iuj m can be completed afterwards. And then the 
Grassman variables can be integrated out. 

To integrate out the Grassman variables, we make a 
change of variables to bring the relevant terms into a di- 
agonal form, with the process similar to the diagonaliza- 
tion of the Hamiltonian of a system of fermions. For this 
purpose, we reexpress the terms containing the Grassman 
variables into a matrix form through an introduction of 
a column matrix of Grassman variables. The resulting 
2x2 matrix is to be diagonalized. Since this 2x2 ma- 
trix is not Hermitian, a proper care must be given to the 
finding of its eigenvalues and corresponding right eigen- 
vectors. With the Grassman variables integrated out, 
the summation over the imaginary frequencies ioj n can 
be performed. 

Lastly, we consider only the real parts of the order 
parameters. For the real order parameters, the imaginary 
part of the form iIm(</j*A — A*tp) in the contribution 
to the action arising from the integration of Grassman 
variables is identically zero. With this imaginary part 
identified to be zero, the imaginary parts of the order 
parameters in the other contributions to the action can 
be integrated out, which provides a reasonable procedure 
for retaining only the real parts of the order parameters. 



3 



Finally, neglecting all the irrelevant prefactors in the 
partition function arising from the above-described al- 
gebraic manipulations, we obtain the following path- 
integral representation of the grand partition function 



dipdA e 



where S is referred to as the action and is given by 

5 = Vln Sinl f A ^-2Vlnf2cosh^ 
^ BAv ^ V 2 

+ 2N8Avip 2 + NBUA 2 

with Av = v — \x and 

TP (i-l | 2 2 TT 2 a2\1/2 

Ek = {i k +g <y5 - U A ) . 
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Here ip and A are the order parameters of zero momen- 
tum and frequency. Notice that the summation over p 
in the first term on the right hand side of Eq. (0) yields 
a factor proportional to the number of molecules with 

We now infer information about the order parameters 
for Bosc-Einstein condensation of molecules and BCS 
condensation of atoms from the above path integral rep- 
resentation of the grand partition function. The integrals 
in the grand partition function can be performed in prin- 
ciple by using the steepest descent method. The station- 
ary point in the order parameter space in this method 
corresponds to the Hartree-Fock-Bogoliubov approxima- 
tion to the order parameters. To obtain an effective 
action, we proceed as follows. With the action on the 
exponential being expanded as a Taylor series of the or- 
der parameters about the stationary point and with the 
terms up to the second order being kept, the integra- 
tion in Eq. Q can be performed. The effective action in 
the Hartree-Fock-Bogoliubov approximation is actually 
obtained by retaining only the stationary-point contri- 
bution to the action. 

Taking partial derivatives of the action with respect to 
ip and A and setting them to zero to find the stationary 
values of ip and A, we obtain the self-consistent equations 
for ip and A 



<P = V 
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(9a) 
(9b) 



The self-consistent equations of the order parameters 
ip and A in Eqs. © are the central results derived in the 
present work and they provide the basis for our follow- 
ing discussions and conclusions. Noticing that U > for 
the Fermi gas of atoms under our current study and the 
value of the wavevector sum in Eq. I|9b|l is positive, we 



immediately realize that there exists only a trivial solu- 
tion for A from Eq. (|9b|> . which implies that fermionic 
atoms can't condense via the mechanism of BCS con- 
densation! However, the self-consistent equation for ip in 
Eq. (|9a|l may have a nontrivial solution. We can thus 
conclude that possible condensation or superfluidity in 
a repulsive Fermi gas of atoms is that of molecules or 
preformed Cooper pairs. 

Algebraically, the nonexistence of a nontrivial solution 
for A is because there is a minus sign in front of A 2 in the 
excitation spectrum of quasiparticles as can be seen from 
Ek = (£fc + g 2 <p 2 - U 2 A 2 ) X I 2 . This minus sign leads to 
a minus sign on the right hand side of the self-consistent 
equation for A. 

Physically, the nonexistence of a nontrivial solution for 
A is because of the absence of an effective attractive two- 
body interaction among atoms. From the Hamiltonian in 
Eq. ||TJ , it is seen that the only possible physical process 
that might generate an attractive two-body interaction 
among atoms is the atom-molecule coupling. However, 
the effect of this coupling is actually pair breaking in that 
a pair is lift out of its condensed state and transformed 
into a molecule. The reverse process does not lead to 
pairing in that it can only be said that a molecule is dis- 
solved into two atoms that occupy two states with oppo- 
site momenta. However, for pairing to occur, there must 
exist an effective attractive interaction between these two 
atoms. 

We now investigate whether an effective attractive in- 
teraction can be produced by the atom-molecule cou- 
pling. In analogy to the phonon-mediated electron- 
electron interaction in metals, it can be explicitly shown 
that the effective interaction between atoms mediated 
through molecules is repulsive. For this purpose, we con- 
sider the following Hamiltonian H = Hq + Hi with 

H = 2Ava*a + ^ £ fe 4 CT c feCT , (10a) 

ha 

Hl = J k E( atcfe T c -fcl + h - c )- (i° b ) 

^ k 

Here for simplicity we included in the Hamiltonian only 
the contributions from zero-momentum molecules and we 
shall consider the effective atom-atom interaction me- 
diated through these molecules. The effective atom- 
atom interaction mediated through nonzero-momentum 
molecules is qualitatively similar. 

To infer the effective atom-atom interaction mediated 
through molecules, we adopt the procedure due to Foldy 
and Wouthuysen |19l \2(\ and make a canonical transfor- 
mation to the above Hamiltonian 



H' = e~ s He s 



with S assumed to be of the form 

S = ^(Aka^ c feT c_ fe + 5 fc acL fe 4 T ) 



(11) 



(12) 
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Expanding H' and keeping terms up to the second order 
in S, we have H' « H + [H, S] + ±[[H , S},S}. The coeffi- 
cients in S are to be determined from the condition that 
Hi + [H 0l S] = 0. The transformed Hamiltonian is then 
approximately given by 

H'hH + ^[Hi,S\. (13) 

The concerned effective atom-atom interaction mediated 
through molecules is to be deduced from the second term 
in the above equation. With only the atom-atom two- 
body terms retained, the effective atom-atom interaction 
is given by 

From the above equation, it is clearly seen that 
the effective atom-atom interaction mediated through 
molecules is repulsive for atoms close to the Fermi surface 
as long as the molecule offset energy Av = v—^iis greater 
than zero. Since atoms close to the Fermi surface play 
a dominant role in BCS condensation and since Av > 
is satisfied in the concerned Fermi gas of atoms ^(J, we 
thus conclude that the effective atom-atom interaction 
mediated through molecules is repulsive. 

For atoms in a repulsive Fermi gas of atoms used in 
the recent experiments m addition to the repulsive 
two-body interaction tunable through the Feshbach res- 
onance, there exists also the above-calculated effective 
interaction mediated through molecules. According to 
BCS theory 0, l^] , the formation and immediate con- 
densation of Cooper pairs of fermionic atoms require that 
the net two-body interaction between atoms be attrac- 
tive. From the fact that both the original two-body in- 
teraction and the effective two-body interaction mediated 
through molecules are repulsive, it follows that the net 
two-body interaction is repulsive, which implies the im- 
possibility of BCS condensation of atoms, in consistency 
with the conclusions drawn previously from the self- 
consistent equations for the order parameters. There- 
fore, the observed condensation on the BEC side can't 
be BCS condensation of atoms. It may be condensa- 
tion of molecules or that of preformed Cooper pairs. No- 
tice that condensation of preformed Cooper pairs dis- 
tingushes from BCS condensation of atoms in that the 
former is a two-step process in which Cooper pairs first 
come into being as bound states of atoms and then con- 
dense, while the latter is a single-step process in which 
the formation and condensation of Cooper pairs occur 
simultaneously. If condensation on the BEC side ob- 
served in the experiment turns out to be that of pre- 
formed Cooper pairs, then the mechanism of the forma- 
tion of these Cooper pairs will certainly shed light on the 
mechanism of high temperature superconductivity. 

To summarize, within the path-integral formalism 
of the grand partition function, we have derived self- 



consistent equations for the order parameters of Bose- 
Einstein condensation of molecules and BCS condensa- 
tion of atoms in a Fermi gas of atoms with a repulsive 
two-body interaction. We have found that there is only 
trivial solution of the order parameter for BCS condensa- 
tion of atoms, from which we conclude that BCS conden- 
sation of atoms is not achievable in such a Fermi gas of 
atoms and that the observed condensation must be that 
of molecules or that of preformed Cooper pairs. The ex- 
plicit calculation of the effective atom-atom two-body in- 
teraction mediated through molecules demonstrates that 
there is no net attractive interaction between atoms and 
therefore confirms our conclusions. 
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